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Q^' 2-particle channel is determined from the energy spectrum obtained in 
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1 Introduction and motivation 

Determination of scattering phase shifts is still considered one of the open 
problems in the nonperturbative lattice approach to quantum field theory. 
Up to now bulk quantities, masses and matrix elements were the central 
objects. Scattering processes involve at least four asymptotic states; in QCD 
each of them is a bound state in terms of the original fields, indeed a very 
demanding situation. However, various contributions to this question exist 
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The problem is particularly pressing, as most of the hadronic states to 
be analysed are resonances and hence not asymptotic states of the field the- 
ory. This makes any rigorous control of finite size effects, which necessarily 
turn up due to numerical limitations, extremely difficult. Furthermore, phase 
shifts provide the only direct comparison to experiment. Resonance param- 
eters are inferred from phase shifts whereas positions of resonance poles in 
the unphysical sheet are plagued by necessary assumptions implicit in the 
phenomenological analysis. 

In a series of papers Liischer [|I[]-[§| introduced a method to obtain scat- 
tering phase shifts from an analysis of the energy spectrum in a finite size 
system. Eventually this method should be used in the full 4-d theory p, 
but before we can hope to attack the determination of phase shifts in lattice 
QCD it is appropriate to deal with simpler models. Liischer and Wolff [Q 
studied the 0(3) non- linear sigma model in d=2, where analytic results for 
the scattering phases exist. That model may be interpreted as describing 
scattering of three mass-degenerate particles with certain symmetries. In or- 
der to discuss a situation of two fields with different masses and potentially 
resonating scattering behaviour we attacked another simple 2-d model [|1^ . 
Here we present a full account of our analysis and results. 

The organisation of this paper is as follows. In the next section we briefly 
review Liischer's approach for the calculation of scattering phases in a two 
dimensional lattice field theory. We also discuss the idea for the numerical 
determination of the scattering spectrum from Monte Carlo data as presented 
in m . Based on this, we formulate a method that allows one to compute the 
contribution of a lattice operator to physical states. In the third section we 
present our model and discuss basic properties and our Monte Carlo simu- 
lation. The fourth section contains results for the scattering parameters. A 
determination of the wavefunction renormalization constants provides bet- 



ter understanding of the observed dependence of the resonance width on the 
bare 3-point couphng. This is followed by a section where we present our 
results for the contributions of the lattice operators to the scattering states. 
We end with some concluding remarks. 

2 Determination of scattering parameters 

2.1 Energy spectrum and phase shift 

Let us start with a brief review of Liischer's approach ll|]-[§|. Throughout 
the paper we consider the dimensional quantities as given in units of the 
appropriate powers of the lattice spacing a. Consider the scattering of two 
bosons of mass m in a. two dimensional system of finite spatial extension 
L, but infinite time extension. In this finite volume the momentum of a 
particle can only have discrete values. Assume that the interaction region is 
localized and smaller than L. Imposing periodic boundary conditions, the 
quantization condition for the relative momentum of two scattering particles 
reads 

2S{kn) + knL = 2mT, neN, (1) 

where 6{k) denotes the phase shift acquired in the interaction region. If its 
functional form is known, one could use this relation to find the quantized 
values kn- On the other hand, given the momentum spectrum, (|^) allows the 
determination of S{k) for each k. This approach may be realized by utilizing 
the dispersion relation for the energy of two particles (of equal mass m) with 
relative momentum k 

W = 2Vm2 + A;2 . (2) 

The energy spectrum in the two particle sector may be computed in a Monte 
Carlo simulation for a lattice with spatial extension L. Indeed, for given L one 
measures a set of discrete values Wn, thus kn and with (|lD the corresponding 
value S{kn). Repeating this procedure for various lattice sizes it is possible 
to find the phase shift for a large set of momentum values. 

The power of this approach lies in the utilization of finite size effects for 
the determination of physical quantities. One will have to respect carefully 
the limitations. The interaction region and the single particle correlation 
length ought to be smaller than the spatial volume, in particular 1/m << 



L. Polarization effects due to virtual particles running around the torus 
should be under control. Lattice artifacts will turn up in O(a^) corrections, 
i.e. for large values of k. For the determination of the energy spectrum 
one should consider correlation functions of a sufficiently large number of 
observables with the correct quantum numbers, capable to represent the 
space of scattering states 0] . 

2.2 Correlation functions 

We compute the connected cross-correlations 

M„^(t) := (iV„*(0)iV„(t))e = (iV„*(0)iV„(t)) - (iV„*)(iV^) (3) 

of operators in the scattering sector. Each Nn,n = 1,2... is defined on a 
single time slice, and t denotes the separation of the time slices. Ideally we 
have to consider a complete set of states. The explicit form of the operators 
will be given in section 4, after introduction of the model. 

The transfer matrix formalism (see e.g. |Tl|) yields the spectral decom- 
position 

oo 

Mnm{t) = j:vil^*vll}e-'^^ . (4) 

1=1 

For simplicity we assume non-degenerate Wi, ordered increasingly. The am- 
plitudes f^'^ = (/|A^„0) are the projections of the states |A^nO) (generated by 
the operators Nn out of the vacuum) on the energy eigenstates (/| of the 
scattering problem. 

For the determination of the energy spectrum from (^ we follow [Q . On 
a finite lattice, in the elastic regime below the four particle threshold (or the 
three particle threshold, depending on the model) there is a finite number of 
energy eigenstates. Only two-particle states contribute and the projections 
onto physical single particle states vanish. Also, in a Monte Carlo calculation 
one can only take into account a finite number r of operators, which, however 
should exceed the considered number of eigenstates. 

One can expect that the spectral decomposition of M{t) is approximated 
rather well by the truncated matrix 

a™(t) = E4n«e-*^' , n,m = l,2...r . (5) 



The generalized eigenvalue problem 

Cit)C^'\t) = \^'\t,to)Cito)C^'\t) , k = l,2,..,r . (6) 

allows one to find the energy eigenvalues Wi efficiently even for not too 
large values of t. In (H) we assume to < t, e.g. to = 1- We require the 
Nn , n = 1,2, ..,r to be linearly independent, thus C(to) is regular and the 
generalized eigenvalue problem is well defined. Indeed, it may be transformed 
to a standard eigenvalue problem for the matrix 

C-i(to)C(t)C-5(to), (7) 

with identical eigenvalues X^'^\t,to) and eigenvectors u'^'^)(t) = C2(to)C (^)- 
The solution for the eigenvalues is 

AW(t,to) = e-(*-*o)'^^ . (8) 

This can be easily verified by inserting it into @, leading to 

B^''\t)C^''\t) = , (9) 

where 

Each row of the matrix B{t) is a linear combination of the r — 1 linearly 
independent vectors \'^''\l = 1,2, ..,k — 1,/c + l,..,r . Thus the space of 
the solutions of {^) has dimension 1, as is needed to construct a normalized 
eigenvector for the corresponding standard eigenvalue problem. The effect of 
the truncation was shown to behave like exp[— (t— to)VFr+i] by a perturbation 
calculation 0. 

To summarize, the determination of the spectrum in the scattering sector 
consists of the following steps. Compute the correlation matrix M{t) for 
various values oft and calculate M~2(to)M(t)M~2 (^q). The dimension r of 
M cannot be chosen too large, because then the inversion turns out to be 
numerically unstable. Finally solve the eigenvalue problem for this matrix, 
and determine the spectrum from the exponential decay (BI). 



2.3 Eigenstates and spectral density 

Like in any Quantum Mechanics application the central problem is to con- 
sider a set of operators with a sufficient overlap with the energy eigenstates. 
We now demonstrate that the eigenvectors q''' of the generalized eigenvalue 
problem contain information about the representation of the physical states 
by the chosen operators Nn- 
Rewriting (|) gives 

^ v«t[(e-*'^' - e-*o^'e-(*-*°)'^^)v« ■ C^'^(t)] = (11) 

for all k. Since the v*^'-'"'' are linearly independent, their coefficients have to 
vanish, yielding v(') ■ C^''\t) = Q \f I ^ k . Thus C^''\t) oc v^*^) and, when 
normalized to unit length, independent of t. Introducing the normalized 
eigenvector Q^' = \^^y \ v^'^) | also the eigenvectors u'^'') = C2(to)C of the 
associated standard eigenvalue problem are time independent. 

We may express the original approximation (^ to M{t) in terms of the 
obtained eigenvectors. 



C^Ut)=ECrC^Mv«pe-*'^' , n,m=l,2...r . (12) 

1=1 

In particular consider the diagonal element 

an(t) = E I cf n v^'M' e-*^' . (13) 

1=1 
For a continuous energy spectrum, 

Cnn{t)= j dWpniWy'"^ , (14) 

and we identify | Cn"* PI v*-'-* poc Pn{Wi), the spectral density of the correlation 
function of operator Nn. We therefore find that | ^^'^ p indicates the relative 
weight of the contribution of operator Nn to the energy eigenstate | I) . 

In the calculation again one has to deal with the correlation matrix M 
and not with the truncated form C. This implies that results obtained from 
^W := M-^(to) u('), where uW are the eigenvectors of M-5(to)M(t)M-^(to), 
are affected by a term 0(e~*'^''+M, where r is the rank of M. 

If energy eigenvalues are degenerate, the procedure does not apply and 
the amplitudes cannot be identified in a unique way. The method then gives 
equal importance for the involved operators. 
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3 Model and Monte Carlo simulation 

3.1 Model 

We choose a model, which describes two hght particles yj that couple to a 
heavier particle t] giving rise to resonating behaviour. For this we couple two 
Ising fields through a 3-point term. The action is given by 

xeA,fM=l,2 

— l^r) 2^ VxVx+fi 
x€A,fi=l,2 

+ ^ J2 Vx^xi^x-fi + ^x+f,) ■ (15) 

a'GA,/i=l,2 

Here, ft denotes the unit vector in direction /i. The values of the fields 
are restricted to {+1,-1}. Quantization corresponds to taking expectation 
values with the measure Z~^ exp {—S) dipdrj where Z denotes the normalizing 
partition function. The sum runs over all sites (a;o, xi) of the euclidean LxT 
lattice A C Zi X Zr and we imply periodic boundary conditions. The 3-point 
term was introduced in a nonlocal but symmetric way, because ip'^ = 1. 

The case g = corresponds to the situation of two independent Ising 
models, each with a 2nd order phase transition at Kcrit = | ln(l + v2) — 
0.44068. In the scaling limit the model describes an interacting boson with 



mass m = — log(tanhK) — 2k |]TT|. Reformulating the Ising model at the 
phase transition as a theory of non-interacting fermions, it was shown that 
the scattering matrix assumes the value —1 independent of the momentum 
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In the coupled case {g > 0), if we identify r] and (p with particle states, 
the term proportional to g gives rise to transitions like rj —>■ ifip rendering 
rj a resonance in the ip channel, when kinematically allowed. The notion of 
mass is then no longer well defined for the rj as it is no longer an asymptotic 
state of the theory. For us, this is just the situation of interest. Henceforth, 
whenever we say "mass" of the r] field we really mean the position where the 
phase shift assumes its resonating value. For given g, "masses" of the fields 
T] and (f may be adjusted by calibrating the hopping parameters n^ and k^. 



3.2 Phase structure 

In the following we give some simple arguments about the phase structure of 
our model in the (k^, k^, g) space (in the thermodynamic limit). 

For ^f = we have the Ising phase transitions for k^ = /tcrit and k^ = Kcrit • 
In the general case the model is even in the ip field. In the symmetric phase 
of the ?7-field there is a symmetry g ^->- —g. 

For g > we first consider the situation with a quenched i] configuration 
and discuss the dynamics of (p. The action then has just a hopping term for 
the yj-field with a locally varying coupling 

uJx,^, = K^- -giVx + Vx+ii) ■ (16) 

For sufficiently small g this is still a ferromagnetic system and the Griffiths 
inequalities apply. Thus all link expectation values are monotonically in- 
creasing with uJx^^- For sufficiently large k^^ there will be an ordered and 
for sufficiently small k^ a disordered phase. This situation (Potts model 



with quenched random bonds) has been recently studied in detail |ll3|, |T4 
where the authors prove that there is a second order phase transition. To 
our knowledge it is unclear, whether this behaviour survives for a dynamical 
Tj field. 

In order to examine the rj dynamics let us discuss two limiting cases. For 
K^ = oo the if spins are strictly ordered. Then the three-point term in (^) 
effectively acts like an external magnetic field of strength g. For non-zero g 
the field t] has long range order and does not undergo a phase transition. On 
the other hand for k^ = 0, y? is completely disordered and the interaction 
term acts like an external random field for the 77-spins. Random magnetic 
fields are expected to lower the effective dimensionality (i.e. raise the critical 
dimension by 2) [|15|-|jl3. Recently it was proved [|1^ that for a quenched 



random magnetic field the field rj remains disordered for all values of the 
hopping parameter. Hence we conclude that there is no phase transition in 
the rj field for finite k^ and Ky, = 0, g > 0. 

We expect therefore that the 1] phase transition extending at g = over 
all values of k^ disappears at g > 0, at least for the edges k^ = 0, cxd. 
Our numerical results for masses and bulk quantities show that the phase 
transition becomes weaker for g > 0, changing into a cross-over like behaviour 
with large, but finite correlation length at the r^-transition. 



To summarize these arguments, one cannot be sure that for (? 7^ the 
model ([T5|) undergoes a second order phase transition for both fields. This 
would indicate that the cutoff cannot be removed completely ior g > and 
that there is no continuum limit possible, in which both fields describe mu- 
tually interacting physical particles with finite mass. We are always dealing 
with an effective model. 

However, we are not presenting a model for a possible continuum field 
theory. All we want is to study the possibility of extracting phase shifts in a 
situation where two particles of different mass have an effective interaction. 
For this aim the lattice is considered more a technical means to approximate 
an underlying continuum model than a method for a rigorous construction 
of a continuum limit. In the range of bare couplings considered the model 
serves our purpose. However, not all combinations of correlation lengths (i.e. 
mass values) can be realized. In particular a study of lattice artifacts is not 
possible for arbitrary small lattice spacing. 

3.3 Cluster algorithm 



For the Monte Carlo simulation of Ising models the cluster algorithm |T8[ is 
a very powerful tool, both for generating new configurations and for defining 



and measuring improved estimators |T9[. In our case with a 3-point term we 



have to modify it. We update the r] and the if field alternately. 

Updating the rj spins: Bonds between neighboured spins of equal 
sign are kept with the probability 1 — exp (— 2/t;^). After identification of the 
connected clusters all clusters are flipped, each with separately determined 
probability 



pfP = l/(l + e-2°(^)) 



"(<^) = o H VxVx{<^x-f, + Vx+il)- (17) 

a:eC,At=l,2 



This method has been suggested |T^ in order to incorporate the effect of 
constant external fields. 

Updating the ip spins: This part of the action is treated like an Ising 
system with the locally varying hopping parameter cjj,^^ introduced in (|T6|) . 
In identifying the clusters we keep bonds between like-sign neighbours alive 
with probability 1 — exp {—2ux,fj,) and flip clusters with probability |. 
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The cluster method is ergodic and fulfills the detailed balance condition; 
for the values of g considered it is a clear improvement over the Metropolis 
updating. 

We study the model at three values of the 3-point coupling g = 0, 0.02 
and 0.04; for later reference we denote these cases by I, II and III. Our 
first step was to tune the hopping parameters. We measured single particle 
propagators to adjust the (y9-mass and find a crude estimate for the energy of 
the 7] resonance. We tried to establish the values m^ ~ 0.19 and ?7i^ ~ 0.5. 
Table 1 gives the couplings where we performed systematic simulations on 
various lattice sizes. 

Throughout this work we use T = 100; the spatial extension L varies 
between 12 and 60. For each set of couplings and lattice size we performed 
typically 2 x 10^ measurements. For the statistical errors we display one 
standard deviation, estimated with the Jackknife method (For a discussion 



see e.g. |^ 



4 Results for the scattering parameters 

For the determination of the energy spectrum a precise knowledge of the 
single particle mass and related finite size effects is important. In this sec- 
tion we first discuss our results for the single particle state, then for the two 
particle channel, in particular the phase shifts. Finally we perform a com- 
putation of wavefunction renormalization constants, to clarify the scaling of 
the resonance width with g. 

4.1 Single particle states 

Let us first study the single particle propagator. The operator of a 93 state 
with momentum 

Pi,^ = 2tiv/L, z/ = -L/2 + 1,...,L/2 (18) 

may be represented through 

1 Y. </'xo,xi exp (ixipi,^) , (19) 
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where A^.^ denotes a tiineslice of A. Its connected correlation function over 
temporal distance t decays exponentially oc exp {—po,u't) defining the single 
particle energy p^y, in particular we have Po^u=o = f^ip- 

The observed mass, as compared to the "real" mass at vanishing lattice 
spacing and infinite volume, incorporates contributions from polarization due 
to self interaction around the torus (decreasing exponentially with L) and 
lattice artifacts due to the finite ultra-violet cutoff (polynomial in the lattice 
constant a). 

In fig. 1 we plot po,i^, obtained for the pure Ising model (I) and L = 50. 
We compare the values with the continuum spectral relation , 



Po,u = V"^' + pIu (20) 

and find deviations 0((ap)^), as expected [Q. We also plot the dispersion 
relation for the lattice propagator of a Gaussian particle with mass m, 

Po,i/ = arcosh(l — cospi,,^ + coshm) , (21) 

and find excellent agreement with the measured values. Thus this particular 
single particle state shows little deviation from the free lattice form. This 
results holds for models (II) and (III) too. 

The Ising model may be considered the limit of a 0^-theory with infinite 
bare 4-point coupling A. We are in the symmetric phase. Thus we expect the 
leading interaction to be due to an effective 4-point term. The finiteness of 
the spatial volume allows self-interaction with one particle running around 
the torus. This leads to a finite volume correction to the particle mass. 
Performing a perturbation calculation for two dimensional ip^ theory along 
the lines of [^ one finds a decrease 



M{L)-m oc AL-3e-'"^ + 0(A^L~5e-2™^) , (22) 

where M{L) is the mass measured on a lattice of spatial size L, and m is the 

infinite volume mass. A 3-point interaction leads to a different exponential. 

Fig. 2, again for model (I), gives our results for ln[A/L(M(L) — m)] 



vs. L and indicates perfect agreement with (22). Indeed the fit to a linear 
behaviour gives a slope of — 1.047(31)?7i for case (I) and — 1.015(25)m for 
(III), in good agreement with the expected value —m. For situation (I) 
this is obvious, because we are in the symmetric phase where there is no 
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3-point coupling. In cases (II) and (III) the rj particle that enters the three- 
point interaction is much heavier than the (p and its contribution is therefore 
suppressed. This mass shift on finite lattices has also been confirmed in 



Monte Carlo simulations of the Ising model in 4 dimensions |22| . 

As already discussed in section 2, the infinite volume mass is explicitly 
known in the Ising model. We compare the measured value at L = 50 for 
g = with this exact result and find a shift of +0.00006(13). We use this 
shift to estimate the infinite volume mass values for the light particles in 
the coupled cases by subtracting it from the values measured at L = 50, 
correspondingly (Anyway this shift lies within the errors). 

Equivalencing the Ising system with an effective field theory of renormal- 
ized fields, the real valued variables $ with the canonical continuum form of 
the kinetic and the mass term for bosons, transform to the Ising spins via 
$ = ^/K(f). The momentum space propagator for the $ fields is 

""^"^ = lufTV' tor p_0 , (23) 

where rrir is the renormalized $ mass and Z the wave function renormaliza- 
tion constant. It may be determined from the susceptibility 

Z = Km'^ Y: (0.,0..)c . (24) 

x1,x2gA 



The renormalized field is given through $,, = ^/y/Z = ^Jk (p/y Z. 

Since the definitions for rj and (p are identical, the symbol stands for 
either spin variable rj or Lp. We use the high temperature results p3[ for the 
susceptibility in the Ising model, to check the correctness of our numerical 
results for Z in the uncoupled case. 

Only Z^p is well defined, since only p is an asymptotic state of the full 
theory. However, in the discussion of our results we want to compare the ob- 
tained value for the resonance width with the tree-level predictions of an ef- 
fective Lagrangian model; for this we need to consider properly renormalized 
effective fields. Only for the sake of this qualitatively interesting comparison 
we determine an "effective" w.f. renormalization constant Z^, aware of the 
ill-definedness of that quantity in case (II) and (III). 
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4.2 Two particle sector 

4.2.1 Operators and phase shifts 

In the 2-particle channel we consider operators with total zero momentum 
and quantum numbers of the 7], 

Ni{xo) = - Y. Vxi,xo , (25) 



1 

Z2 



^.(^o) = 7^ E e^"^(--^^)^.„,,^,„,, , (26) 



xi,yieAxQ 



277 (j - 2) 
with pj = , J = 2, 3, . . . 

and measure all cross-correlations @. The operators Nj^i describe two 
(y9-particles in the CM system with relative momentum pj. Due to the in- 
teraction they do not correspond to eigenstates of our model. Indeed they 
are eigenstates of the Gaussian model of free bosons. However, if the set is 
complete, the diagonalization of section 2 provides the necessary information 
on the energy spectrum of this channel. 

The number of operators considered in (|^) should be chosen larger than 
the number of states in the elastic regime 2?n^ < W < Arriy, and not larger 
than L/2 to be linearly independent. A larger set provides a better rep- 
resentation of the eigenstates but enhances the numerical noise. We work 
with between 4 and 6 operators depending on L (Only as a testing case we 
once consider 8 operators for a particular run at L = 60, g = 0.02). For the 
product M~2 (to)M(t)M~2 (to)- we choose to = 1; for our statistics, usually 
the diagonalization becomes unstable around t = tmax — 5 — 8, depending 
on L. The values Wa are determined by averaging log[Aa(t)/AQ(t + 1)] for 
1 < t < tmax- For each lattice size we determine the low lying part of the 
energy spectrum. Each value corresponds, via the spectral condition (^, to 
a value of k and with (|I|) to a phase shift 6{k). 

For g = we have the situation of two independent Ising systems. As 
discussed above, the critical Ising model in d=2 has an S-matrix equal to 
— 1 ||l2|, and we therefore expect a phase shift 6ismg = ~f (mod vr) in the 



scaling regime. In fig. 3 we exhibit the energy spectrum (a) for a system 
of non-interacting Gaussian fields (curves due to (|l]) with phase shift 5 = 0) 
and (b) as obtained in our simulation of the Ising model. The curves for 
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the two particle states in fig. 3b are due to (^ ioi 6 = — | and we find 
excellent agreement. The single particle state t] (horizontal line at W = 
0.5) is completely decoupled from the (fip states and we therefore observe 
degeneracies in the energy levels at certain values of L. Fig. 4 gives the 
results for the phase shift S{k) determined from the data for the Ising model 
(fig. 3b) as a function of the dimensionless momentum k/m^p. 

Already here we identify two sources for the relative magnitude of the 
error bars. Small values of k result from values of W near 2'm^p, at large L. 
Close to this threshold, due to the functional dependence of @j the statistical 
error in W transforms into a relatively larger error of k and thus 5{k). Higher 
A;, on the other hand, stem from the large values of the energy, where the 
statistical fluctuations are intrinsically larger. 

Figs. 5a and 5b exhibit results for g ^ Q. For most of the points the errors 
are smaller than the symbols, typically 0.5% of the energy value. The degen- 
eracies in the energy spectrum have now disappeared, resulting in avoided 
level crossing with a gap that grows with g. The corresponding phase shift 
data are given in fig. 6. 

One should keep in mind that points at neighboured values of k may come 
from quite different values of L and different branches of energy levels. In 
order to emphasize this situation we have chosen corresponding symbols to 
denote energy values on one branch, i.e. for one value of n. Each branch, 
from small to large L contributes to the whole range of k values; it starts at 
larger values of W (or k) with a phase shift close to 7r/2, passes through the 
resonance value at the plateau and approaches — 7r/2 towards smaller k. The 
overall consistent behaviour is impressive. 

The high energy values, corresponding to large k may feel the lattice cut- 
off effects 0{a?) as well as possible mis-representation of the energy eigen- 
states by the considered operators. We discuss the second problem in a 
subsequent section. Within the accuracy of our data none of these effects 
seems to be a problem. 

4.2.2 Resonance parameters 

The observed overall behaviour has a simple interpretation. A resonance may 
be parametrized with the effective range approximation 

k 
— cot (5Res = a-hk"^ . (27) 
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At the resonance S^es — f + r ^2 ' °^ 



m« = 2v/(m^)2 + (a/6), Yr = ^^- . (28) 

Adding this resonance phase to the Ising model background phase the spec- 
trum in the resonating channel should approximately obey 

(^totaK^;) = Rising + (^Res ( A;) = -arctan——, (29) 

I 2 2r?7r 

W = 2y(m^)2 + (--5totai(A;) + -^)' , n = 0,l,... .(30) 

Here n denotes the various energy levels in relation (|l|) . We superimpose this 
theoretically expected behaviour of the spectrum on the measured points in 
figs. 3-6. 

The parameters rriji = rrij^ and T^ = F^ may be determined from a direct 
fit to the shoulders in the energy curves. However, this fit turns out to be 
not very stable and has much larger statistical errors than a determination 
from the resulting phase shifts as discussed below. We fit the values of tan S 
in the region —1.45 < 6 < 1.45 to a straight line in k"^ according to (|29|) . 
This gives the effective range parameters a, b and our estimates for m^, F^ 
in table 1. The resulting x^ indicates that the statistical errors of our phase 
shift data and the errors of W may be underestimated (cf. [0). 

The ratio of the resonance widths for (HI) and (II) is smaller than the 
value 4, the square of the ratio of the couplings g that one naively expects 
from the Born approximation. In doing that, however, we should first write 
down an effective Lagrangian with properly renormalized effective fields rjr 
and ifr, with kinetic parts and a three-point term. For sake of such a com- 
parison we have to determine the wave function renormalization constants 
Z^ and Zrj, although the second is ill-defined as mentioned earlier in sec. 4.1. 
From this information we may estimate an effective three-point coupling 

9cS = 9 ^ Z^JZr, , (31) 

K Ik V 

which may be a more suitable number to enter the tree-level formula. 

With this naive procedure we obtain Z^, Z^ and ^^efr for all three values 
of the bare coupling g (see table 1). Both Z^ and Z^ are smaller than 1, even 
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in the uncoupled case. These values decrease for increasing g. The values 
for ^eff are 0.0568(6) and 0.0938(8) for g = 0.02 and g = 0.04 respectively. 

With these values, the ratio (5feff(0.04) / 5fcff(0.02)) = 2.73(11) is closer 
to the observed ratio r^(0.04) / r,,(0.02) = 2.13(26). The naive argument 
based on tree-level perturbation theory appears to be at least qualitatively 
correct. However, for a better identification of our model with an effective 
field theory, one has to take into account further interaction terms of higher 
order in the fields. In fact the preferred definition of a renormalized 3-point 
coupling is just from the observed (physically well defined) width. 

5 Representation of the physical states 

Is the number of lattice operators considered for the determination of the 
energy levels in the scattering channel sufficiently large? To analyse the 
representation of the physical states by the considered operators ( p5|) we 
compute the vectors Q' = M~2(to) u'^'\ normalized to unit length. Then 
[Cn^j provides a measure on how much iV„ contributes to the physical eigen- 
state I /). As discussed in sect. 2.3 there should be no dependence on t. This 
is consistent with the numerical results and therefore we average the values 
for all available t. 

Fig. 7 shows the relative weights (Cn^) for / = 1, 2, 3 and n = 1, 2, 3 for 

case (III). Due to the normalization of Q ' the individual contributions can 
assume values between and 1. If the three displayed operatois are sufficient 
to represent the considered state, their contributions should add up to 1. 

Since these values are displayed for all considered lattice sizes, a remark 
concerning the number r of components in Q', which is equal the rank of 
M, is appropriate. For small L, the lattice is statistically "nervous", and the 
inversion of M is possible only for a small number r of operators (e.g. r = 4 at 
L = 12, while r = 8 at L = 60). At small L, for the lowest energy eigenstate 
|1) in the scattering sector mainly A^^i contributes (see the circles in fig. 7a). 
Only Q is considerably different from zero. Hence for the lowest state |1) 
it makes no difference whether the vector Q ' is normalized at length 3 or 4, 
because the 4th entry is very small anyway. For the higher states |/) more 
than four entries of Q ' may be necessary. However, due to the numerical 
instability of the inversion at small L, only four entries are available. Thus 
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the results for the representation of the higher physical states appear shifted 
for small L. For larger L this technical problem vanishes: The results for at 
least the three lowest eigenstates are reliable. 

Compare the energy spectrum of fig. 5b with fig. 7. Let us first discuss 
the lowest energy state |1) , represented by circles in fig. 5b, with respect 
to fig. 7a. For small L, Wi is close to the resonance energy. There the 
state is completely dominated by the operator A^i (circles in fig. 7a), which 
describes an r^-particle at rest. For increasing L, Wi decreases and approaches 
the 2-particle threshold. This manifests itself in a drastic decline of the 
contribution of A^^i, accompanied by an increase of the amplitudes for two 
(/9-particles at rest (operator A^25 triangles) and with relative unit-momentum 
(operator N-^, squares). For lattice size larger than 20, the energy Wi is 
already too far below the resonance energy so that the energetically higher 
7] at rest cannot contribute much to |1). The state is built out of N2 and N^ 
(with energies as shown in fig. 3a) almost exclusively. 

For the second lowest state |2) (fig. 7b) the situation is reversed. As seen 
in fig. 5b, for small L the energy W2 is high above the resonance energy, so 
the ry-state Ni (circles) does not contribute. The state is dominated by the 
lowest lying 2-ip operators. For increasing L, W2 decreases and as it crosses 
the resonance energy, the contribution of A^^i peaks. Leaving the resonance- 
plateau, state 1 2) is again dominated by the energetically favourable operator 
iVa. The contribution of A^4 is substantial as well, but not shown in the 
picture. 

Fig. 7c gives the contributions to |3). In the considered L- region it 
does not show any dramatic behaviour. For L < 50 the energy is too high 
for Ni and N2 contributions and the state is built out of A^3 and higher 
operators. Only for the largest L the ?7-amplitude grows. For small L the 
remark concerning the rank of the correlation matrix made earlier applies. 
Between L = 12 and L = 16 there is a discontinuity in the N3 contributions, 
due to the discussed truncation effect. The diagonalization at L = 16 includes 
one operator more than that at L = 12. This operator has nonvanishing 
contributions and therefore affects the A^s amplitude via the normalization 
ofC^^'). 

In the case (I) of two uncoupled Ising models we get a similar picture. 
Ni then measures the particle mass in the completely disconnected r] system 
and does not affect the (p states and vice versa. The lowest eigenstates states 
of fig. 3b turn out to be mixtures of the ipip-staies Nn , n > 2. 
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In particular, it is easy to identify misrepresentation of energy eigenstates 
in terms of the considered operators. At higher energies one observes a shift 

of the weight factors [Cn^) towards higher operators A^„. Fig. 8 gives an 
example of a situation, where the center of weight has moved to relatively 
high values of n. It is obvious, that more operators should be included for a 
good representation of that state. 

Numerically we observe, that a truncation of the correlation matrix M{t) 
leads to a rise of the higher energy levels. However, in the data presented 
here and used in the determination of the phase shift this effect is under 
control. For the low lying energies it is negligible and for the largest energy 
level it is still smaller than the statistical error. 

6 Concluding remarks 

In our Monte Carlo study we have established, that Liischer's suggestion for 
determining phase shifts is indeed a very reliable method, at least in d=2. 
Our simple model describes a system of two particles with a three-point 
coupling. The coupling constants have been chosen such that one of the 
states may become a resonance. Indeed we find resonating phase shifts in 
this channel. As a byproduct we also obtain the phase shift of a pure Ising 
model. 

The following remark concerns just the particular model studied here: As 
discussed, for g ^ we cannot expect to perform a continuum limit. Most 
likely it is not possible to increase the correlation length arbitrarily keeping 
the ratio m,j/m^ fixed. Numerically we observed, that the ?7-contribution to 
the peak in the specific heat becomes less pronounced for increasing values 
of g (This is the technical reason for not studying even broader resonance 
widths). Therefore, in this model, it might not be possible to study lattice 
artifacts in a more systematic way. 

The study of the representation of scattering states by the operators en- 
tering the correlation functions sheds light on the problem of constructing 
good approximations for the energy eigenstates. In particular it allows to 
control unwanted truncation effects. 

For small L one is confined to the first plateau only. Translating a L- 
independent constant energy W = m^ into a phase shift we get a variety 



of 5- values, all for a constant k = \J'm'^ — Arn^, just what we expect for 
a zero width resonance: a stable particle state. For a naive estimator of 
the resonance energy one therefore might consider just the single particle 
operator, ignoring the states of higher energy coming from above. However, 
increasing the lattice size would not improve the result (In fact, the resonance 
is no asymptotic state of the theory). There can be no avoided level crossing, 
as there are no 2-particle states. Thus there will be no information on the 
resonance width. On the other hand, ignoring the r^-state and considering 
only ipif states, is insufficient, too. It is of crucial importance to include all 
operators (with the quantum numbers of that sector) with energies close to 
the considered energy eigenstates. For further discussions, in particular with 
regard to the situation in QCD see 0. 

How can we understand the limit (? ^ 0? The energy gap becomes 
smaller and the level crossing region narrower. The phase shift approaches 
step function behaviour. In the Ising limit the phase shift is then compatible 
with a constant (mod vr). 
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Tables 

Table 1: This table summarizes the parameters for our simulations together 
with the measured values of masses and resonance width (in units of the 
lattice constant). We also show wave function renormalization constants and 
the effective 3-point couplings. 



Model 


I 


II 


III 


9 


0. 


0.02 


0.04 


Knp 


0.3948 


0.3897 


0.3700 


K^i 


0.3268 


0.3323 


0.3700 


m^ 


0.189921 


0.18965(13) 


0.19014(13) 


rrirf 


0.499545 


0.5008(4) 


0.5112(3) 


r. 


0. 


0.0047(4) 


0.0100(3) 


z. 


0.722(2) 


0.706(4) 


0.671(4) 


Zfj 


0.849(3) 


0.816(5) 


0.619(4) 


9cs 


0. 


0.0568(6) 


0.0938(8) 
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Figures 



Fig. 1: The dispersion relation of our single particle propagator (circles, 
model (I) for L = 50) approaches the continuum behaviour (broken curve) for 
small pi^i^; it follows nicely the lattice propagator of a free Gaussian particle 
(full curve). 

Fig. 2: We plot our values for log[(M(L) — m)\/T\ vs. L, model (I); 
these results for the single particle mass follow the expected exponential 
dependence on the spatial size (^2]). 



Fig. 3: (a) The 2-particle energy levels for a free Gaussian model {S = 0) as 
functions of the lattice size L. The dotted line indicates the energy of a stable 
rj with rrirf = 0.5. (b) The measured energy levels for the ipip and r] channel 
for the coupling constant values g = vs. the spatial lattice extension L. 
The full curve shows the theoretical expectations for the values from (0) for 
Rising = ~f ) the horizontal line denotes the r] mass level. The dashed lines 
indicate the 2- and the 4-particle thresholds. 

Fig. 4: The results of the phase shift for model (I), i.e. corresponding to 
the energy data in fig. 3b, vs. k/m^. The full curve gives the theoretical 
expectation ^ismg = ^f • "^^^ dashed line indicates the 4-particle threshold. 

Fig. 5: The measured energy levels for the (fip channel for the coupling 
constant values (a) g = 0.02 and (b) g = 0.04, vs. the spatial lattice extension 
L. The full curves show the theoretical expectations for the values from (|^) 
and (5totai as discussed in the text, eg. (POD . The dashed lines indicate the 2- 
and the 4-particle thresholds. 

Fig. 6: The phase shift as determined from the energy levels in figs. 5a and 
5b vs. k/m^p] the full curves denote the theoretical expectations as discussed 
in the text. The dashed line indicates the 4-particle threshold. 

Fig. 7: Contributions of the operators Ni,N2 and N^ (circles, triangles, 
squares) to the energy eigenstates (a) |1) , (b) |2) and (c) |3), for case (III), 
displayed as a function of the lattice size L. 
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Fig. 8: Contributions Kn^M of tlie considered operators A^„ to tlie energy 
eigenstate |6) at L = 60 for model (II). 
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